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Abstract 

We study the hadron-quark phase transition in neutron star matter and the structural properties of liybrid stars 
using an equation of state (EOS) for the quark phase derived with the Field Correlator Method (FCM). We make 
use of measured neutron star masses, and particularly the mass of PSR J1614-2230, to constrain the values of the 
gluon condensate G2 which is one of the EOS parameter within the FCM. We find that the values of G2 extracted 
from the mass measurement of PSR J1614-2230 are fully consistent with the values of the same quantity derived, 
within the FCM, from recent lattice QCD calculations of the deconfinement transition temperature at zero baryon 
chemical potential. The FCM thus provides a powerful tool to link numerical calculations of QCD on a space-time 
lattice with neutron stars physics. 
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1. Introduction 

Neutron stars, the compact remnants of su- 
pernova explosions, are unique natural labora- 
tories to explore the phase diagram of quan- 
tum chromodynamics (QCD) in the low tem- 
perature T and high baryon chemical potential 

region [1,2]. In this regime nonperturbative 
aspects of QCD are expected to play a crucial 
role, and a transition to a phase with deconfined 
quarks and gluons is expected to occur and to 
influence a number of interesting astrophysical 
phenomena [3-9]. 

Recent high precision numerical calculations 
of QCD on a space-time lattice at /i^ = (i.e. 
zero baryon density) have shown that at high 
temperature and for physical values of the 
quark masses, the transition to quark gluon 
plasma is a crossover [10] rather than a real 
phase transition. 

Unfortunately, present lattice QCD calcula- 
tions at finite baryon chemical potential are un- 
realizable by all presently known lattice meth- 
ods (see e.g. [11]). Thus, to explore the QCD 



phase diagram at low T and high /i6, it is nec- 
essary to invoke some approximations in QCD 
or to apply some QCD effective model. 

Along these lines, different models of the 
equation of state (EOS) of quark matter, as 
the bag model [12] or the Nambu Jona-Lasinio 
(NJL) model [13,14], have been intensively used 
by many authors to calculate the structure of 
strange stars [15-17], or the structure of the 
so called hybrid stars, i.e. neutron stars with 
a quark matter core. These EOS models are 
expected to be reasonable at very large density, 
but they crumbles in the density region where 
quarks clusterize to form hadrons, i.e. in the 
region where the deconfinement phase transi- 
tion takes place. In addition, the bag model 
and the NJL model, as other QCD effective 
models, can not make predictions in the high T 
and zero /it. region, and thus can not be tested 
using present lattice QCD calculations. 

Recently the deconfinement phase transi- 
tion has been described using an EOS of quark 
gluon plasma derived within the Field Corre- 
lator Method (FCM) [18,19] extended to finite 
baryon chemical potential [20-22]. The FCM 
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is a nonperturbative approach to QCD which 
includes from first principles the dynamics of 
confinement. The model is parameterized in 
terms of the gluon condensate G2 and the large 
distance static quark-antiquark (QQ) poten- 
tial Vi. These two quantities control the EOS 
of the denconfined phase at fixed quark masses 
and temperature. The main constructive char- 
acteristic of the FCM is the possibility to de- 
scribe the whole QCD phase diagram as it can 
span from high temperature and low baryon 
chemical potential, to low T and high /ib limit 
[20-22]. 

A very interesting featiirc of the FCM is that 
the value of the gluon condensate can be ob- 
tained from lattice QCD calculations of the de- 
confinement transition temperature Tc, at zero 
baryon chemical potential. Thus we have an ef- 
ficacious tool to directly link lattice QCD sim- 
ulations and neutron star physics. 

To explore this link is the main purpose of the 
present work. In particular, we will investigate 
the possibility for the occurrence of the quark 
deconfinement transition in neutron stars and 
the possibility to have stable hybrid star con- 
figurations using the FCM for the quark phase 
EOS and a relativistic mean field model [23] for 
the EOS of the hadronic phase. 

2. EOS ofthe Quark Phase 

The quark matter equation of state we used 
in this work is based on the Field Correlator 
Method [18,19]. Recently this method has been 
extended to the case of nonzero baryon density 
[20-22] making possible its application to neu- 
tron star matter. 

The main advantage of the FCM is a natural 
explanation and treatment of the dynamics of 
confinement in terms of Color Electric D^(x), 
Df{x) and Color Magnetic D"{x), D^{x) 
gaussian correlators [19]. 

contributes to the standard string tension 
cr^ through [20]: 

a^ = \j D^{x) <fx. (1) 

The string tension vanishes as goes to 
zero at T > Tc [20], and this leads to deconfine- 
ment. The correlators have been calculated on 
the lattice [24] and also analytically [25]. 

Within the FCM the quark pressure Pq, for 
a single flavour, read [20] : 



^il^ =:^\.^''y — y — ) + M J' — )](2) 

where 

f 1 

4>M)= / du - — , — rrTT'^^^ 

J V w + 2^ (exp vw + z^^ — a + 1) 


u = niq/T and Vi is the large distance static 
QQ potential: 

l/T 00 

V,{T) = J dT{l - tT) J dxX^f (VX^TT^) .(4) 


The nonperturbative contribution to Df (x) is 
parametrised as [19]: 

Df{x)=Df{0)exp{-\x\/X) (5) 

where A is the vacuum correlation length. Fol- 
lowing Ref. [20], we use the value A = 0.34 fm 
which has been determined in lattice QCD cal- 
culations [26]. 

In this formalism Vi in Eq.(4) is indepen- 
dent on the chemical potential (and so on the 
density) . This feature is partially supported by 
lattice simulations at small chemical potential 
[20,27]. In the present work, the value of Vi at 
T = has been considered as a model parame- 
ter. 

The gluon contribution to the pressure is [22]: 

00 

= ^ / dxx' , \y,, . • (6) 

37r^ y exp(x-h|^)-l 

In summary the total pressure of the quark 
phase is given by: 

^s^^.+ EA-^G^. (7) 

The last term in Eq.(7) represents the vac- 
uum energy difference between the quark 
and hadronic phases, in the case of three fla- 
vor {u, d, s) quark matter [20], and G2 is the 
gluon condensate. The latter quantity has 
been determined, with large uncertainty, using 
QCD sum rules [28] to be in the range G2 = 
(0.012 ± 0.006) GeV^. In the present work, the 
value of G2 has been considered as a model 
parameter. We used the following values of the 
current-quark masses: rriu = rrid = 5 MeV and 
TOg = 150 MeV. In summary, the quark matter 
EOS has two parameters: G2 and Vi =Vi{T = 
0). 
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3. Neutron star strucure 

In this section we show the results of our 
calculations of hybrid stars structure. To this 
purpose we integrate the well known Tolman, 
Oppenheimer and Volkov relativistic hydro- 
static equilibrium equations (see e.g. [29]) to 
get various stellar properties for a fixed EOS. 
For the hadronic phase we consider /3-stable 
nuclear matter, and we make use of a nonlinear 
relativistic mean field model in the parametri- 
sation GMl given in [23]. The GMl model 
can be considered a representative realistic 
nuclear EOS in the sense that it fits the em- 
pirical saturation properties of nuclear matter, 
does not violate causality at high density, and 
is compatible with present measured neutron 
star masses. All the results presented in the 
following have been obtained using the Gibbs 
construction [30] to model the hadron-quark 
phase transition. 

In Fig. 1 we report the stellar gravitational 
mass M (in unit of the solar mass Mq = 
1.99 X 10"^^g) versus the central baryon num- 
ber density pc (left panel) and the mass versus 
radius R (right panel) in the case of pure nuclc- 
onic stars (continuous line) and of hybrid stars 
for different and taking Vi = 0.01 GeV. 
We obtain stable hybrid star configurations 
for all the considered values of the gluon con- 
densate, with maximum masses ranging from 
M^ax = 1.44 M© (case with G2 = 0.006 GeV^) 
to Mmax = 2.05 Mq (G2 = 0.016 GeV'^). No- 
tice that the hybrid star branch of the stellar 
equilibrium configurations shrinks as G2 is 
increased. This behaviour is different with re- 
spect to the one found in Ref. [31] where the 
stability window of hybrid star configurations 
was restricted between 0.006 GeV^ < G2 < 
0.007 GeV*. 

The properties of the maximum mass configu- 
ration for hybrid star sequences varying G2 are 
summarized in table 1. 

In Fig. 2 we plot the quark-hadron phase 
transition boundaries in /3-stable nuclear mat- 
ter as a function of G2 and with Vi = 0.01 GeV. 
The onset of the deconfinement transition (i.e. 
the onset of the quark-hadron mixed phase) oc- 
curs at the baryon number density pi , and the 
pure quark phase begins at p2. Also shown is the 
central baryon number density p^^'' of the max- 
imum mass hybrid star (dotted-dashed line). 
Stable hybrid star configurations have central 
densities pc located in the region of the P-G2 



plane between the lower continuous line and the 

dotted-dashed line, i.e. pi < Pc < P^^^ ■ Notice 
that p^y^ > p2 when the gluon condensate is in 
the range 0.006 GeV^ < G2 < 0.0077 GeV. For 
these G2 values all hybrid stars with a central 
density in the range P2 < Pc ^ Pc^^ possess a 
pure quark matter core. Finally the horizontal 
dashed line represents the value of the central 
baryon number density p^^ of the maximum 
mass pure nucleonic star. 

In Fig. 3 wc draw the maxinnuxi mass M„i.„x 
for hybrid stars (continuous line) and the mass 
Ml = M{pi) (dashed line) of the star with cen- 
tral baryon number density pi corresponding to 
the onset of the mixed phase. These two quan- 
tities are plotted as a function of the gluon con- 
densate G2 and taking Vi — 0.01 GeV. Stable 
hybrid star configurations correspond to the re- 
gion of the M-G2 plane between the continu- 
ous and the dashed line. Stellar configurations 
in the region below the dashed line Mi do not 
possess any deconfined quark matter in their 
center (pure nucleonic stars). 

To compare our results with measured 
neutron star masses, we report in the same 
Fig. 3 the values of the masses of the fol- 
lowing pulsars: PSR B1913-hl6 with M = 
1. 4398 ± 0.0002 Mq [32]; PSR J1903-F0327 with 
M = 1.667 ± 0.021 Mq [33]; and PSR J1614- 
2230 with M = 1.97 ± 0.04 M© [34]. 

The mass of PSR J1614-2230 gives the 
strongest constraint on the possible value of 
the gluon condensate. In fact, using the lower 
bound of the measured mass of PSR J1614- 
2230, wc get G2 > 0.0129 GeV*. Thus for 
values of the gluon condensate in the range 
0.0129 GeV < G2 < ~ 0.018 GeV*, 
PSR J1614-2230 is a hybrid star, whereas 
PSR B1913-hl6 and PSR J1903-F0327 are pure 
nucleonic stars. In the above specified range for 
the gluon condensate, G2 is defined by the con- 
dition Mi(G2) = 2.01 Mq, the upper bound of 
the measured mass of PSR J1614-2230. Thus 
for G2 > G*2 PSR J1614-2230 is a pure nucle- 
onic star. 

To explore the influence of the large distance 
static QQ potential on the stellar properties, 
we have considered a quark phase EOS with 
Vi = 0.10 GeV. Once again we get stable hybrid 
star configurations for all the considered val- 
ues of G2 , with maximum masses ranging from 
M^ax = 2.00 Mq (case with G2 = 0.006 GeV*) 
to Mmax = 2.25 Mq (G2 = 0.0016 GcV*). 
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Thus an increase of the value of Vi reduces the 
extension of the hybrid star branch, shifts it 
to larger densities and produces hybrid stars 
with a larger maximum mass. In this case we 
found that the calculated Mmax is compatible 
with the lower bound of the measured mass of 
PSR J1614-2230 for all the considered values 
of the gluon condensate {i.e G2 > 0.006 GeV^). 

We also considered stellar models with Vi = 
0.50 GcV and Vi = 0.85 GeV. In these two 
cases no phase transition occurs in neutron 
stars [i.e p\ > p^^ the central density of the 
maximum mass pure nuclconic star), thus in 
this case PSR J1614-2230 would be a pure 
nucleonic star. 

4. Lattice QCD calculations and 
measured neutron star masses 

Within the FCM the dcconfinement transi- 
tion temperature Tc at /ib = reads [20] 



(8) 



with ao = (37r^/768)^/^ in the case of three fla- 
vors. 

In their analysis the author of Ref. [20] assume 

Vi{Tc) = 0.5 GeV, thus in Eq. (8) is a simple 
function of G2, and is represented in Fig. 4 by 
the curve labeled Vi (Tc) = 0.5 GeV. This result 
can hence be compared with lattice QCD calcu- 
lations of Tc giving the possibility to extract the 
range of values for the gluon condensate com- 
patible with lattice results. This comparison has 
been done by the authors of Ref. [20] , and it is 
done in the present work in Fig. 4, where we 
consider recent lattice QCD calculations of Tc 
[36,35]. Details to the specific lattice QCD cal- 
culations are given in the Fig. 4 caption. As one 
can see, the comparison with lattice QCD calcu- 
lations of Tc restricts the gluon condensate in a 
rather narrow range G2 = 0.0025-0.0050 GeV*. 

Next to verify if these values of G2 are com- 
patible with those extracted in Sec. 3 from 
hybrid star calculations and measured neutron 
star masses, we need to relate the parameter 
Vi = V\{Q), entering in the zero temperature 
EOS of the quark phase, with 14 (Tc) in Eq.(8). 
To this end, one can integrate Eq.(4) using the 
non perturbative contribution (5) to the color 
electric correlator Df{x) and assuming that 
the normalization factor ff (0) does not de- 



pend on temperature. The latter assumption is 

supported, up to temperatures very near to Tc, 
by lattice calculations [24]. Therefore one gets: 

»,,,r, = ;.,(0){i-|f + l(i + 3f 

Thus Vi(Tc) = 0.5 GeV corresponds to Fi(0) = 
0.85 GeV to be used in the T = EOS of the 
quark phase. In this case, as we found in Sec. 3, 
no phase transition occurs in in neutron stars 
{i.e pi > p^^) for all the considered values of 
G2 ■ Thus for these values of the EOS parameters 
PSR J1614-2230 would be a pure nucleonic star. 

We can also evaluate the FCM transition 
temperature at fii, = corresponding to the 
case Vi{0) = 0.01 GeV used in Sec. 3 for hy- 
brid star calculations with the T = FCM 
equation of state. To this purpose we solve 
numerically Eqs. (8), (9) and we obtain the 
results represented in Fig. 4 by the curve la- 
beled Vi = 0.01 GeV. The comparison of these 
results with lattice QCD calculations [36,35] of 
Tc restricts the gluon condensate in the range 
G2 = 0.0103 0.0180 GeV*. Coming now to the 
astrophysical constraints on the gluon conden- 
sate, the vertical grey line in Fig. 4 represents 
the lower limit for G2 which is compatible, 
in the case Vi(0) = 0.01 GeV, with the lower 
bound of the measured mass of PSR J1614- 
2230 (see Sect. 3). 

A similar analysis can be done for the case 
Vi{0) = 0.10 GeV. Now the comparison be- 
tween the FCM transition temperature at pb = 
(curve labeled Vi = 0.10 GeV in Fig. 4) and 
lattice QCD calculations of the same quantity 
gives G2 = 0.0085-0.0153 GeV'', whereas one 
gets G2 > 0.006 GeV^ from the lower bound of 
the measured mass of PSR J1614-2230. 



5. Conclusions 

In this letter we have studied the hadron- 

quark dcconfinement transition in /3-stable nu- 
clear matter and the structural properties of 
hybrid stars using an EOS for the quark phase 
derived from the Field Correlator Method ex- 
tended to finite baryon chemical potential. We 
obtained stable hybrid star configurations for 
all the values of the gluon condensate fulfilling 
the condition pi{G2) < p^^{G2), i-e. the de- 
confinement transition can occur in pure nucle- 
onic stars. 



.(9) 
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We have established that the values of the 
gluon condensate extracted within the FCM 
from lattice QCD calculations of the decon- 
finement transition temperature at /Xb = are 
fully consistent with the value of the same 
quantity derived by the mass measurement of 
PSR J1614-2230. The FCM thus provides a 
powerful tool to link numerical calculations of 
QCD on a space-time lattice with neutron stars 
physics. 
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Fig. 1. Stellar gravitatonal mass M versus central 
baryon number density pc (left panel) and versus stel- 
lar radius R (right panel) for hybrid stars for several 
values of the gluon condensate G2 (reported in GeV* 
units) and for Vi = 0.01 GeV. The continuous line in 
both panels refers to the pure nuclconic stars, i.e. com- 
pact stars with no quark matter content. 
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Fig. 3. Gravitational maximum mass for hybrid stars 
(continuous line) and gravitational mass Mi (dashed 
line) of the star with central baryon number density 
pi corresponding to the onset of mixed quark-hadron 
phase as a function of the gluon condensate G2 and for 
Vi = 0.01 GeV. 
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Fig. 2. Quark-hadron phase transition boundaries in 
/3-stable nuclear matter as a function of the gluon 
condensate G2 and for V\ = 0.01 GeV. The onset of 
quark-hadron mixed phase occurs at the baryon num- 
ber density pi, and the pure quark phase begins at 
P2. Also shown is the central baryon number density 
p^^^ of the maximum mass hybrid star. The horizon- 
tal dashed line represents the value p^^ of the central 
baryon number density of the maximum mass pure nu- 
cleonic star. 
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Fig. 4. (Color online) Deconfinement transition tem- 
perature Tc at pny = 0. The curve labeled with 
Vi(Tc) = 0.5 GeV reproduces the FCM results of 
Ref. [20] for a fixed value Vi{Tc) = 0.5 GeV of the 
large distance static QQ potential. The curve labeled 
with Vi = 0.01 GeV (Vi = 0.10 GeV) corresponds 
to the transition temperature at /if, = obtained 
solving numerically Eqs. (8) and (9) for the case 
Vi(0) = 0.01 GeV (Vi(0) = 0.10 GeV). The horizon- 
tal heavy and thin lines represent respectively the cen- 
tral value and the error esimate of lattice QCD calcu- 
lations. In particular, the (red) continuous lines refer 
to the calculations [35] of the HotQCD collaboration 
Tc = (154 ± 9) MeV; the (blu) short-dashed lines re- 
fer to the calculations [36] of the Wuppertal-Budapest 
collaboration Tc = (147 ± 5) MeV. Finally, the vertical 
green line represents the lower limit for G2 which is 
compatible with the lower bound of the measured mass 
of PSR J1614-2230 for the case Vi(0) = 0.01 GeV. 
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G2 [GeV*] 


Mmax/MQ 


[fm-3] 


R [km] 


0.006 


1.44 


1.55 


9.54 


0.012 


1.89 


0.77 


12.55 


0.016 


2.05 


0.75 


12.66 



Table 1 

Properties of the iiiaximuin mass configuration for hy- 
brid stars as a function of the gluon condensate G2. 
The results are relative to the case Vi = 0.01 GeV. The 
maximum mass configuration for the pure nucleonic 
star sequence is: Mmax = 2.33 M©, p^^ = 0.87 fm~* 
and R = 11.70 km. 
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